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ABSTRACT 

The current status of lattice calculation of weak matrix elements for heavy 
quark systems is reviewed. After an assessment of systematic errors in 
\ present simulations, results for the B meson decay constant, the B parame- 

ON ' ter Bb and semi-leptonic heavy-to-light and heavy-to-heavy transitions are 

discussed. The final topic are lattice results for heavy baryon spectroscopy. 
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1. Introduction: Lattice Approach to Heavy Quark Systems 



Heavy quark systems play an important role in the study of the less well-known el- 
ements of the CKM matrix, which are central to our understanding of the origin of CP 
violation, and may also contain information of physics beyond the Standard Model. Set 
against this background, the lattice formulation of QCD provides a framework for the 
calculation of hadron masses and weak matrix elements from first principles. Since the 



lattice approach is intrinsically non-perturbative, one may tackle the large theoretical 



uncertainties due to the strong interaction in weak decay amplitudes. 

Lattice QCD replaces space-time by a four- dimensional, euclidean, hypercubic lattice 



4h 

of size Lr • T. The sites are separated by the lattice spacing a, which acts as an UV cut-off. 
One problem encountered in current simulations is that typical values of a' 1 lie in the 
range 2-3.5 GeV. Therefore one expects that discretisation errors ("lattice artefacts") will 
distort the results already for charm physics. Also, b quarks cannot be studied directly, 
since their mass is above the UV cut-off. 
^ ■ Several methods are being used to circumvent this problem. It is now customary for 

simulations in heavy quark physics to cancel the leading discretisation error by employing 
so-called improved actions and operators 0, or by absorbing it into a rescaling of the quark 
fields i. Quantities can then be computed around m charm and extrapolated to the b quark 
mass. Alternatively, one can use the static approximation and perform the simulation at 
infinitely heavy quark mass, using the leading term of the expansion of the heavy quark 
propagator in I/tuq. Finally, one can employ a non-relativistic formulation of the QCD 
Lagrangian. It is obvious that none of these methods is entirely satisfactory, but that 
they provide complementary information about heavy quark systems. 

Apart from lattice artefacts, the main systematic errors in lattice simulations in- 
clude the effects due to neglecting internal quark loops by using the so-called quenched 
approximation. The normalisation of lattice operators is another source of systematic 
uncertainties; due to the explicit breaking of chiral symmetry by the regularisation proce- 
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dure, lattice operators are in general related to their continuum counterparts via (finite) 
normalisation constants, whose numerical values are usually not known very precisely. 
Also, the explicit breaking of chiral symmetry and Lorentz invariance leads to mixing 
with higher dimension operators. Hence, the connection between the matrix element of a 
continuum operator O and its lattice counterparts is in general given by 

0)^ = J2 Z i(9 2 )(Oi) lm + 0(a), (1) 

i 

where the Zj's are the appropriate normalisation and matching factors. 

Finally, lattice estimates of dimensionful quantities are subject to uncertainties in the 
lattice scale. They arise from the fact that different quantities which are used to set 
the scale a -1 [GeV] give different results. This is closely related to using the quenched 
approximation, since loop effects are not expected to be the same for different quantities. 

2. Leptonic B Decays and B° — B° Mixing 

The decay constants of heavy-light pseudoscalar and vector mesons, fp and fy, are 
related to the matrix elements of the lattice axial and vector currents trough 

(0|4 att (0)|P) ~M P fp/Z A , (0\v; att (0)\V)~ e] M*(fvZv)-\ (2) 

where Za, Zy are the normalisation constants of the currents. The matrix elements as 
well as the pseudoscalar and vector masses Mp, My are extracted from mesonic two- 
point correlation functions through a fitting procedure. HQET predicts scaling laws for 
the combination fp\/~Mp, which is expected to behave like a constant as Mp — > oo. 
Furthermore, in the infinite mass limit, the HQ-spin symmetry predicts that pseudoscalar 
and vector decay constants become degenerate (up to short-distance corrections). Hence 

U(M) = Ml I {l + £a s (M)} = 1 + 0(1/M) (3) 

It has been shown that U (M) indeed approaches one in the infinite mass limit i. However, 
lattice studies have revealed that the 1/M corrections to the scaling laws are large; for 
the decay constant they amount to about 15% at the B mass and about 40% at the mass 
of the D meson. 

In Table 1 we show the results for the decay constants fo, /d s and fp from various 
lattice calculations using relativistic heavy quarks.^ Expressing the results for the deca 



•ay 



constant at the individual values of a in terms of a common hadronic scale ro — 0.5 fm 
one can extrapolate the dimensionless quantity fpr^ to the continuum limit, a — > 0. 
Repeating the procedure using lattice data for f w , one obtains the continuum result from 

f P [MeV] = 132 MeV {f P r,/f n r ) j^. (4) 
The results from this analysis using relativistic heavy quarks read 

f D = 205 ± 45 MeV, f Ds = 220 ± 50 MeV 

f b = 170 MeV, f B Jf B = 1-13 ±0.14. (5) 
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Table 1: Lattice estimates for the pseudoscalar decay constants from different collaborations. 
All data were obtained in the quenched approximation using relativistic heavy quarks. 



Collab. 


a [fm] 


Jd [MeV] /D a [MeV] / B [MeV] f B Jf B 


MILCS 





182(3)(9)(22) 198(5)(10)(19) 151 (5) (16) (26) 1.11(2)(4)(8) 


RLS 5 

1JUJ 1 — 1 


n 059 


?08(9V35Vl 2"l 230(5Vl OVl 9"l 1 87f1 0V34V1 1 1 Kf,) 


T A 1MT N 
Li/VIM Li u 


U.Uoo 




229^/J _ 16 /!oU(4J _ 22 
186(29) 218(15) 


ukqcdI 

n 


0.069 
0.083 


185 +3 + 4 7 212(4) + 4 7 160(6) til 1-22 +3 
199(15) tH 225(15) 1™ 176(25) ±f| 1.17(12) 


ELC0 


0.051 


210(40) 230(50) 205(40) 


jlqcdI 


0.059 
0.078 


216(17) 240(17) 182(16) 
206(12) 237(14) 192(11) 


PCW0 

pcwi 




0.083 


170(30) 180(50) 1.09(2)(5) 
198(17) 209(18) 



The B meson decay constant together with the B parameter B B is of great importance 
for the study of B° — B° mixing. The renormalisation group invariant B parameter B B 
is defined via B B = a^y 2 '^ (5° | L (ji) | S°)/|/|M|, where L (fi) is the AB = 2 
four fermion operator. In a recent study, estimates for the B parameter using the static 
approximation were obtained 

B Bd = 1.02 tl tl B Bs = 1.04 tt tl (6) 

where the first error is statistical, and the second is an estimate of systematic errors. 
The authors attribute a further systematic error of 15-20% to the uncertainty in the 
perturbative matching factors. The above result indicates that SU(3)-flavour breaking 
effects are small for B B , whereas they can be quite sizeable in the case of f B . 
These results can be applied to the ratio of B° — B° mixing parameters x s /xd 



x 1 = r_ B± VB ± M B± f^B^ pi = \V 
x d r Bd fi Bd M Bd p Bd B Bd \V td \ 2 K >\V t 



2 

(7) 
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In conjunction with the experimental value Xd = 0.71(6) li^l, the above results can also 
be used to predict x s , provided l^-j is constrained using global fits0. Choosing Bk = 



0.8 ± 0.2, and taking our lattice estimate f B = 170 ± 55 MeV one obtains 

x s = 13.4 ±4.0 t W 3 3 7 tl% (8) 

where the first error mainly reflects the error in the ratio f B Jf B , the second is due to the 
uncertainty in the actual value of f B , and the third arises from the uncertainty in Bk- 
Clearly, much more precise values of f B are needed in order to predict x s more reliably. 



tWe use a convention in which f„ = 132 MeV. 
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Table 2: Lattice results for form factors relevant for semi-leptonic B — ► tt£u£, ptui decays. 
All results are obtained using propagating heavy quarks with the leading discretisation errors 
subtracted, except for ELC. 



CoUab. 


a [fm] 


/+(0) V(0) A^O) A 2 (0) 


ELC0 


0.051 


0.30(14)(5) 0.37(11) 0.22(5) 0.49(21)(5) 


APE0 


0.083 


0.35(8) 0.53(31) 0.24(12) 0.27(80) 


ukqcdS§ 


0.069 


0.23(2) 0.27 tl ±\ 0.28 tl t\ 


GSsi 


0.059 


0.50(14) t\ 0.61(23) tl 0.16(4) ±™ 0.72(35) ± w 7 



3. Semi-leptonic B decays: heavy-to-light and heavy-to-heavy transitions 

Recently, there has been much activity in studying the decays B — > iriVe, plug, which 
can be used to extract V u b- The matrix elements for these decays are parametrised in 
terms of form factors, e.g. for B — > n£T>£ 

(n\V,\B) = (p B +P,),U(q 2 ) - (PB-P,),f-(q 2 ), (9) 

where q = pb — Pn is the momentum transfer. For the decay B —>■ ptu^ one has additional 
form factors V, Ai, A 2 and A 3 . An important ingredient in the analysis of these decays 
is the observation that for infinite heavy quark mass, HQET predicts scaling laws for the 
form factors near maximum momentum transfer g max , i.e. at leading order in 1/M@ 

/ + (<?LJ ~ ^ 1/2 > n4)~^. Aiifijj-M-W... (10) 

Due to the limitations imposed by the lattice spacing, simulations are not yet suited 
for a direct computation of form factors for B — > itiv^ piv^. One rather obtains the 
form factors for typical lattice momenta \p\ < 1.5GeV/c and for heavy quark masses 
in the region of charm. Hence, the "generic" semi-leptonic heavy-to-light transition in 
current lattice simulations is D — > KiVe, and typical momentum transfers lie in the 
range — 0.8GeV 2 /c 2 < q 2 < 1.7GeV 2 /c 2 . Lattice results for form factors relevant for 
semi-leptonic D decays can be found in a recent review 0. 

In order to make predictions for B — > tc£Vi, pip£, one needs to extrapolate the lattice 
form factors to the mass of the b quark using the above scaling laws. Since the range of 
accessible lattice momenta is rather restricted, i.e. \p\ <C m&, one obtains the form factors 
in a narrow range of q 2 near g max . Therefore, in order to determine the g 2 -behaviour of 
form factors for B — > 7r£p£, pive, or their values at q 2 = 0, one cannot avoid introducing 
a certain model dependence in the lattice results: assuming vector pole dominance is not 
reliable, since the accessible range of q 2 is rather narrow, and lattice data can presently 
not distinguish between different types of pole behaviour. In an alternative procedure, 
one first interpolates the lattice data to q 2 = for quark masses around m charm . But in 
order to extrapolate the resulting form factors at q 2 to the mass of the b quark, one needs 
to guess its leading scaling behaviour in the heavy mass, which, in contrast to eq. (|10D, 
cannot be obtained from HQET. Several methods have been applied, and Table 2 lists 
the results for form factors at q 2 = from various groups. It has been noted that the 
model dependence could be avoided in the framework of light-cone sum rules0, where one 
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Table 3: Lattice results for the slope parameter p 2 of the Isgur-Wise function, the applied 
parametrisation and, where quoted, estimates for \V c b\ using measured decay rates. In the last 
column, we also list the form factor on which the estimate is based, or which method was used 
to formulate heavy quarks. 



Collab. 


a [fm] 


o 2 j o 2 Par. VwJ 

ru,d rs 1 u1 ' 1 ' CD\ 




BSS@ 


0.059 


1.21(26)(33) lin. 0.044(5)(7) 


h+ 


UKQCD@ 


0.069 


0.9 ±1 ti 1.2 tl t\ BSW 0.037(1)(2) t\ 


h+ 


UKQCD0 


0.069 


1.1(5) 1.2 +3 BSW 0.037(3)(5) 




LANL@ 


0.083 


0.97(6) BSW 


K 


MO0 


0.170 


0.95 quad. 


static 


Ken EH 


0.083 


0.41(2) 


static 


HM§ 


0.083 


0.70(17) lin. 


NRQCD 



obtains a leading scaling behaviour of F(0) ~ M~ 3 / 2 for all form factors. This argument 
has so far not been directly applied in lattice simulations. 

It has also been suggested^ to use lattice form factors to calculate the differential 
decay rate for the exclusive decay B — > p + £~T> e beyond the region of charm production. 
One can then avoid the difficult determination of F(q 2 = 0), and a model independent 
extraction of \V u b\ is possible using experimental data for the differential decay rate dT/dq 2 

\Vu b \- 2 ^- 2 OC {\H+(q 2 )\ 2 + \H~{q 2 )\ 2 + \H°(q 2 )\ 2 } = A* (l + B(q> - &J) , (11) 

where if ± (g 2 ) and H°(q 2 ) are combinations of the form factors Ai(q 2 ), A 2 (q 2 ) and V(q 2 ). 
The term A 2 (l + B(q 2 — q^ a&x )) parametrises long-distance hadronic dynamics. In fact, A 2 
plays the same role as the Isgur-Wise function in the case of heavy-to-heavy transitions. 
Here, lattice data can be used to determine the normalisation, whereas in the heavy- 
to-heavy case, the overall normalisation is provided by the HQ-symmetry. In a recent 
study0, the authors obtain A 2 = 21 ± 3 GeV 2 and B = (-8 tj) 10~ 2 GeV" 2 . 

Heavy-to-heavy transitions like B — > DlVt, D*iVe are also parametrised in terms of 
six form factors, h+,h-,hv,hA 1 ,hA 2 an d h,A 3 , which are functions of u>, the product of 
4-velocities of the B and D mesons. HQ-symmetry relates the six form factors to one 
universal form factor, £(u;), called the Isgur-Wise function, which is normalised at zero 
recoil, £(1) = 1. The form factor h + {uj) is related to via 

h + (u) = (1 + /3+(u) + 7+ (w)) (12) 

where /3+(u) parametrises radiative corrections between HQET and full QCD, and 7+(cu) 
denotes (unknown) corrections in the inverse heavy quark mass. 

On the lattice one typically obtains the form factors from a ratio of the relevant matrix 
elements at u and at zero recoil (u — 1). Thereby, some of the systematic effects cancel, 
most notably the normalisation of the axial and vector currents. In order to obtain an 
estimate for £(a>), the known radiative corrections can also be subtracted. Lattice data 
for the form factors h + and h Al have been used to test the HQ-symmetry H 2 !. It has 
been found that h + (u) shows only a weak dependence on the heavy quark mass, and that 
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Table 4: Masses and mass splittings of heavy baryons from various collaborations. 





Collab. 


h = charm 
Latt. [GeV] Exp. [GeV] 


/i =beauty 
Latt. [GeV] Exp. [GeV] 




ukqcdN 
pcwil 


2.27 +3 tl 2.285(1) 


5.64 tl tl 5.641(50) 
5.728±0.144±0.018 


n h 


UKQCDE3 

ukqcdN 
ukqcdS 

UKQCD0 


2.46 +3 +5 2.453(1) 
2.44 + _\ t\ 2.530(7) 
2.41 t\ t\ 2.468(4) 
2.68 t\ t\ 2.704(20) 


5.77 „q +4 5.814(60) 

5.78 +5 6 + 3 5.870(60) 
5.76 tl tt 

5.99 tl 




Collab. 


Latt. [MeV] Exp. [MeV] 


Latt. [MeV] Exp. [MeV] 


Ah-P 


UKQCD0 
ELci 

UKQCD0 
PCW§ 


408 t£ tit 417(1) 
564 ± 88 ± 18 


359 tH 362(50) 
720 ±160 + 13 ° 

Aor . +100 +30 
- 90 -30 

458 ± 144 ± 18 


- A h 


UKQbDt-3 


inn +50 +13 icn/o\ 
190 - 43 -13 1 69 ( 2 ) 


1 KT +52 +11 1 70/1 1 \ 

157 _ 64 I u 173(11) 


n h — n /i 
O* — O' 


UKQCD0 


-17 tH -2 77(6) 

-20 tH -3 83 
90+6+3 

zo -14 -2 


-6 tit t\ 56(16) 

7 +4 +1 
' -8 -1 

+2 +1 

~ 8 -5 -1 



the Isgur-Wise function extracted from both h + and hA 1 is compatible within statistical 
errors. This may be interpreted as a manifestation of the HQ-spin-flavour symmetry. 

In order to extract \V c b\ from the experimentally measured decay rate for B — ► D*£T7£ 
using lattice data, one first needs to parametrise One particular parametrisation is 

&sw(«) = ^ exp {-(2p 2 - l)£=i} , (13) 

where p 2 is the slope of at zero recoil. Other parametrisations, based on linear, 
quadratic and pole ansdtze, have also been studied. Lattice data have so far not revealed 
any significant dependence on the chosen parametrisationii'ii. 

Table 3 contains lattice results for the slope, either for massless spectator quarks (p 2 d ) , 
or in the case that the spectator quark is a strange quark (p 2 ). Also, the parametrisation 
and the estimates of \V c b\ are given. Lattice results for the slope are consistent among 
the different collaborations, although the errors are still large. Recently, preliminary 
results for semi-leptonic — > A c decays have been shown, which will be valuable in 
further studies of the HQ-symmetry. 



4. Heavy Baryon Spectroscopy 

There has been increased experimental and theoretical activity in the study of baryons 
containing one heavy quark. Lattice calculations can make predictions for the masses of 
heavy baryons, which are extracted from the exponential fall-off of correlation functions of 
suitably chosen interpolating operators. The results from two recent studies B ^ are listed 
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in Table 4, and one observes good agreement with experimental data. Mass splittings 
involving heavy baryons have also been studied by a number of groups, and the results 
are also given in Table 4. The agreement between experimental and lattice data for the 
A— Pseudoscalar and S — A splittings is rather good. However, the spin splittings, in 
particular for S* — S c , are definitely inconsistent with experiment. This is attributed to 
a convolution of lattice artefacts and the quenched approximation. 
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